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ИССЛЕДОВАНИЕ ПАРАМЕТРОВ ВОЛНЫ ПРОРЫВА 
ПРИ РАЗРУШЕНИИ РЕЗЕРВУАРОВ ОБЪЁМОМ ДО 30000 М? 
В ЛАБОРАТОРНЫХ УСЛОВИЯХ 


Показана возможность использования лабораторно-экспериментального метода исследования для решения ком- 
плекса вопросов, связанных с нахождением оптимальных геометрических параметров различных по конфигурации ограж- 
дений для предупреждения распространения пожара пролива при разрушении типовых резервуаров номинальным объ- 


ёмом до 30000 mê. 


Ключевые слова: лабораторное моделирование, резервуар, разрушение, волна прорыва, критерии подобия. 


работе [1] представлен разработанный 

специалистами кафедры пожарной безо- 
пасности технологических процессов Академии ГПС 
МЧС России лабораторный стенд для исследования 
параметров волны прорыва, образующейся при ква- 
зимгновенном разрушении вертикального стального 
резервуара (РВС), а также воздействия волны про- 
рыва на ограждения различного конструктивного 
исполнения, в том числе на защитную стенку резер- 
вуаров типа «стакан в стакане» (РВСЗС). При этом 
основное внимание было уделено описанию условий 
гидродинамического подобия исследуемого процес- 
са, конструкции лабораторного стенда и обработке 
полученных экспериментальных данных с оценкой 
чисел Фруда и Рейнольдса на соответствующих 
участках движения потока, показавших их удовлет- 
ворительную сходимость с данными, полученными 
при обработке материалов натурного эксперимента 
по разрушению РВС-700 м? с водой [2]. 

При решении задач с использованием лабо- 
раторно-экспериментального метода исследова- 
ний, связанных с нахождением оптимальных геоме- 
трических параметров различных по конфигурации 
ограждений для РВС и РВСЗС и их устойчивостью 
к воздействию волны прорыва, требуется модели- 
рование потока жидкости при квазимгновенном 
разрушении наиболее широко распространённых 
на объектах нефтегазовой отрасли резервуаров 
номинальным объёмом до 30000 më. С учётом 
поставленной цели на разработанном лабора- 
торном стенде (рис. 1) были проведены дополни- 
тельные опыты с оценкой параметров потока при 
разрушении модельного резервуара, выполненного 
в масштабе 1:130 (предельный для данной конструк- 
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ции стенда масштаб) по отношению к натурному 
РВС-30 000 m°. 

Исходный модельный резервуар имеет диа- 
метр d, = 0,35 M и высоту h, = 0,30 м и жёстко за- 
креплён на направляющих стенда, поэтому при моде- 
лировании натурного резервуара типа РВС-30000 m° 
с сохранением базового диаметра (d, = 0,35 м) 
соблюдение условия геометрического подобия обес- 
печивали соответствующим уровнем жидкости в мо- 
делируемом резервуаре. 

В соответствии с требованиями ГОСТ 31385- 
2016 «Резервуары вертикальные цилиндрические 
стальные для нефти и нефтепродуктов. Общие тех- 
нические условия» рекомендуемые геометрические 
параметры для типового РВС-30 000 m? составляют: 
внутренний диаметр 45,6 м; высота стенки 18,0 м. 
С учётом максимального заполнения резервуара 
на 95 % уровень жидкости в нём составит 17,1 м. 
Тогда в выбранном масштабе 1:130 уровень жидко- 
сти в модельном резервуаре составит 0,132 м. 


/ yji 


Рисунок 1. Общий вид лабораторного стенда 
с модельным резервуаром и защитной стенкой 
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Поскольку в исследуемых процессах, связанных 
с образованием и движением волны прорыва при раз- 
рушении резервуара, главную роль играют силы тяже- 
сти (критерий Фруда, Fr), то должно обеспечиваться 
соблюдение равенства критерия Фруда в соответству- 
ющих точках потоков по трассе растекания [3—6]: 


2 2 
Fr = U, = Uy d (1) 
BH, EuP 


где U „ и, — характерная скорость волны прорыва при 
разрушении резервуара, м/с; H „ h ,— характерная вы- 
сота волны прорыва при разрушении резервуара, м; 
5, = 5, – ускорение свободного падения, м/с? (здесь 
и далее индексы «н» и «м» показывают отношение ве- 
личины соответственно к натуре и модели). 

Другим необходимым условием моделиро- 
вания является соблюдение турбулентности пото- 
ков как на натурном объекте, так и в модели, то есть 
Ке, > Re, (здесь Ке, — критерий Рейнольдса в моде- 
ли; Re, — критическое значение критерия Рейнольдса 
при моделировании открытых потоков, принимае- 
мое, как правило, не менее 10“ [6, 7]): 


Re, = “is, (2) 


re V, — кинематическая вязкость жидкости (в данном 
случае для воды при 20 °С - 1,01:10% м?/с [8]). 

Вследствие того что процесс движения волны 
прорыва при разрушении резервуара представляет 
собой быстроизменяющееся неустановившееся дви- 
жение жидкости, то необходимо соблюдение также 
тождественности критерия Струхаля (критерия го- 
мохронности) [3—6]: 


sh = Vele Ш (3) 
H, h, 
где Т,, Ё, — характерное время процесса, с. 


Помимо указанных выше условий, влияние по- 
верхностного натяжения должно быть настолько ма- 
лым, чтобы оно не препятствовало образованию волн, 
в связи с чем скорость потока со свободной поверхно- 
стью на модели должна быть более 0,23 м/с [3]. 

В целях проверки выполнения рассмотрен- 
ных выше условий гидродинамического подобия на- 
ряду с лабораторными исследованиями, описание 
которых приводится ниже, было выполнено числен- 
ное моделирование процесса разрушения натурного 
РВС-30 000 м? с водой в программном комплексе для 
инженерных расчётов LSDYNA. Подробное описание 
модели, в основе которой лежит численное решение 
систем уравнений гидродинамики и уравнений де- 
формируемого твёрдого тела в связанной простран- 
ственной постановке, представлено, в частности, 
в работах [9, 10]. На рисунке 2 приведены характер- 
ные фрагменты течения жидкости при моделировании 
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Рисунок 2. Характерные фрагменты течения жидкости 
при моделировании процесса разрушения PBC-30 000 m? 


процесса разрушения рассматриваемого типа ре- 
3epByapa, при этом фиксация основных параметров 
потока (Н,, Ц.) в определённый момент времени от 
начала разрушения резервуара (Т.) проводилась co- 
ответственно на расстояниях L = 10, 15, 20, 25 и 30 м 
по трассе растекания от стенки резервуара на постоян- 
ной высоте в 1,3 м от нулевой отметки (табл. 1). 


Таблица 1 


Результаты численного моделирования параметров 
потока при разрушении РВС-30 000 m° 
и оценки критериев подобия 


Sh, 


8,71 32,38 


Ce [er ас 
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T 
=} 
о 


Первоначальный уровень 
воды в резервуаре (0,132 м) 


Вертикальная стенка 
с участками [, 


Фрагмент модельного 
резервуара 


Рисунок 3. Общий вид участка стенда перед 
проведением опытов по оценке высоты потока 
в соответствующие моменты времени на расстояниях /, 


Для сравнительного анализа критериев подо- 
бия была выполнена серия лабораторных исследова- 
ний параметров потока при разрушении модельного 
резервуара в масштабе 1:130. Методика проведения 
экспериментов первой серии опытов по оценке высо- 
ты потока заключалась в следующем. На основании 
лабораторного стенда, выполненного из листового 
алюминия толщиной 0,01 м, непосредственно от стен- 
ки модельного резервуара по трассе растекания по- 
тока обустраивалась вертикальная стенка, конструк- 
тивно выполненная из листового алюминия высотой 
0,3 ми длиной 0,5 м. Стенка имела толщину 0,0005 м, 
что не препятствовало при проведении опытов про- 
хождению вдоль неё потока. Для фиксации высоты 
потока в соответствующие моменты времени на стен- 
ку предварительно была нанесена миллиметровая 
разметка с выделенными участками на расстояниях 
1, = 0,077; 0,115; 0,154; 0,192 и 0,231 м, соответству- 
ющих в выбранном масштабе моделирования рассто- 
яниям L ‚а также первоначальным уровням жидкости 
в модельном резервуаре, равным 0,132 м (рис. 3). 

Далее с помощью кнопки плавного спуска 
(см. рис. 1) модель резервуара с пневматическим 
поршнем, давление в котором создавалось воз- 
душным компрессором, прижимали к основанию 
лабораторного стенда с уплотнительным резино- 
вым кольцом, обеспечивая герметичность системы. 
Резервуар заполняли водой до требуемого уров- 
ня, предварительно установив внутри его обечайки 
вкладыш для имитации квазимгновенного разруше- 
ния внутреннего резервуара. Вкладыш представляет 
собой свёрнутый в цилиндр и вставленный внутрь 
обечайки лист поликарбоната толщиной 0,0006 м 
и высотой 0,132 м, который при заполнении резер- 
вуара водой плотно прижимается к его обечайке 
по всей образующей. При нажатии кнопки резкого 
подъёма (см. рис. 1) с помощью пневматического 
поршня происходило быстрое поднятие обечайки 
модельного резервуара. Освобождённый вкладыш 
под действием напора воды свободно раскрывался 
на 180°, имитируя разрушение резервуара по вер- 
тикальному шву. Образовавшийся поток жидкости 
свободно растекался по основанию лабораторного 
стенда (рис. 4). 


14 


Пожары и чрезвычайные ситуации: предотвращение, ликвидация 


Процесс образования и распространения пото- 
ка по трассе растекания регистрировался цифровой 
фотокамерой марки Nikon 1]1, создающей 5-секунд- 
ные видеоролики с замедленной съёмкой формата 
HVGA 640х240/400 кадров в секунду. Полная рас- 
кадровка видеороликов позволяла устанавливать 
начальный кадр разрушения модельного резервуара 
с образованием волны прорыва, а также фиксацию 
количества кадров при прохождении потоком соот- 
ветствующих участков по трассе растекания l „ что в 
конечном счёте давало возможность оценивать высо- 
ту волны л, в любой момент времени ё, изучаемого 
процесса. 

Важно отметить, что время от момента раз- 
рушения модельного резервуара до снижения уров- 
ня жидкости на основании лабораторного стенда 
в зоне возможной его инструментальной оценки 
(более 0,005 м) не превышало 0,5 с. 


t, = 0,228 c (T, = 2,6 с) 


t, = 0,491 с (T, = 5,6 с) 


Рисунок 4. Характерные фрагменты течения жидкости 
при разрушении модельного резервуара 


Поэтому для оценки местной скорости потока 
использование приборного оборудования, например, 
микровертушки или трубки Пито с дифференциаль- 
ным манометром, практически не представляется 
возможным. Данные обстоятельства обусловлены, 
во-первых, особенностями изучаемого процесса: бы- 
строизменяющееся неустановившееся движение с 
растеканием на открытой поверхности, в данном слу- 
чае по основанию лабораторного стенда. Во-вторых, 
конструктивными особенностями этих приборов. 
Так, использование гидрометрической микровертуш- 
ки позволяет определять только среднюю скорость 
в потоке за определённый промежуток времени, 
как правило, превышающий 1 с. Время же отклика 
цифровых дифференциальных манометров, исполь- 
зуемых в современных модификациях трубок Пито 
для измерений местной скорости, превышает 0,5 с. 
Таким образом, в настоящих опытах для измерения 
местной скорости на соответствующих участках 
трассы растекания потока использовалась специаль- 
но изготовленная динамическая трубка Пито, пред- 
ставляющая собой Г-образную трубку, конструктив- 
но выполненную из органического стекла наружным 
диаметром 0,007 м и толщиной стенки 0,001 м. Вер- 
тикальная часть трубки имела высоту 0,3 м, а участок 
трубки, загнутый под углом 90°, имел длину 0,07 м. 

Вычисление скорости в месте установки труб- 
ки производилось по формуле [3] 


и, = K [2gH,,, (4) 


где К — коэффициент, зависящий от вязкости жидко- 
сти и конструктивных особенностей трубки и опре- 
деляемый при тарировании трубки; ie — скоростной 
напор, определяемый по уровню жидкости в трубке 
от поверхности потока, м. 
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В целях нахождения коэффициента К была 
выполнена тарировка изготовленной трубки посред- 
ством сравнения результатов измерений средней ско- 
рости потока, выполненных соответственно трубкой 
и имеющейся в наличии тарированной гидрометри- 
ческой микровертушкой цифровой модернизирован- 
ной (ГМЦМ) типа «МИКРО-01». На рисунке 5 показан 
общий вид изготовленной трубки Пито и гидрометри- 
ческой микровертушки. 

Измерение скорости жидкости проводилось на 
участке трубопровода длиной 1 м и диаметром 0,05 м. 
С одной стороны трубопровода была установлена за- 
глушка, а в верхней его части смонтирован шаровой 
Г-образный кран, подключённый к системе централь- 
ного водоснабжения. Посредством плавного откры- 
вания крана в трубопровод подавалась вода, создавая 
поток глубиной hy = 0,010; 0,015; 0,020 u 0,025 m co- 
ответственно. 

Свободно вытекающая с другой стороны тру- 
бопровода вода попадала в приёмный лоток и да- 
лее в систему канализации. В верхней центральной 
части корпуса трубопровода заранее было выпол- 
нено окошко размером 0,1х0,03 м, куда поочерёдно 
в центр создаваемого потока устанавливались микро- 
вертушка и трубка. 

Принцип действия датчика ГМЦМ основан на 
электролитическом способе формирования импуль- 
сов. Между двумя электродами, погружёнными в 
водный поток, вращается лопастной винт, изменяя 
электрическое сопротивление цепи. Импульсы напря- 
жения подаются в регистрирующее устройство для 
обработки. Регистрирующее устройство состояло из 
блока формирования импульсов, блока задания ко- 
эффициентов градуировочного уравнения, тактового 
генератора, блока управления и вычисления, блока 
счёта и дешифрации, блока индикации. Регистри- 
рующее устройство было специально адаптировано 


ГМЦМ 


МИКРО 


СКОРОСТЬ М/С. 


6 


Рисунок 5. Общий вид трубки Пито (a) и гидрометрической микровертушки (6) 
для измерения скорости потока жидкости в лабораторных условиях 
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Таблица 2 


Результаты измерений средней скорости потока в трубе микровертушкой и трубкой 


Измерения ГМЦМ 


h ‚м 


пр” 
М, об. Ирмцм м/с, Н 


по формуле (5) 


м 


Измерения трубкой 


,M/c**, 


р по а (4) 


„м/с, Comal um 


по ри (4) 


o ИСИ ПОЗ ПОСЛ ПОС o СИ 
ее ыы Е 


* Без учёта К; ** с учётом К = 0,972. 


к условиям проводимых экспериментов. Принимая 
во внимание номер используемого лопастного винта 
(№ 84), а также градуировочный штекер, соответ- 
ствующий времени осреднения (2 с), использовалась 
следующая формула для определения средней скоро- 
сти потока: 


илим = 0,03787 + 0,0166, (5) 


где n = N/AT, N – показания цифрового индикатора 
регистрирующего прибора (число импульсов); АТ = 
= 2,0 с — время осреднения. 

В таблице 2 представлены результаты измере- 
ний средней скорости потока в трубе соответственно 
микровертушкой и трубкой, а также вычисленный та- 
рировочный коэффициент К. 

Методика проведения экспериментов во вто- 
рой серии опытов по оценке местной скорости потока 
тарированной трубкой в целом аналогична методике 
оценки высоты потока, описанной выше. При этом 


|, = 0,154 м; 
t, = 0,412 с; 
н, = 0,079 м 


вместо стенки по трассе растекания потока на соот- 
ветствующих расстояниях /, от стенки модельного 
резервуара при неизменной высоте в 0,01 м (1,3 м 
для натурного резервуара) устанавливалась трубка. 
Процесс также снимался цифровой фотокамерой. 
Полная раскадровка видеороликов позволяла уста- 
навливать начальный кадр разрушения модельного 
резервуара с образованием волны прорыва, а так- 
же фиксацию в соответствующий момент времени 
значения скоростного напора на соответствующих 
участках потока. На рисунке 6 представлены ха- 
рактерные фрагменты измерений трубкой скорост- 
ного напора на соответствующих участках по трас- 
се растекания потока при разрушении модельного 
резервуара. 

Далее, используя формулу (4), находили со- 
ответствующую местную скорость потока. В таб- 
лице 3 представлены результаты обработки экспе- 
риментальных данных, а также численные значения 
соответствующих критериев подобия. 


Рисунок 6. Характерные фрагменты измерений трубкой скоростного напора 
на соответствующих участках по трассе растекания потока 
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Таблица 3 
Результаты экспериментального определения 
параметров потока при разрушении 
модельного резервуара и оценки критериев подобия 


И, м 


м? 


и» М/С Fr Sh 


M M 
i 
> 
D 


Таблица 4 


Сравнительный анализ критериев 
гидродинамического подобия потоков 


Численные значения Экспериментальные данные 


10,70 


18,18 0,272 17,89 | 23393 
0,115 
23,78 0,360 23,99 | 19928 


29646 


32,38 0,412 31,80 | 19046 
0,154 


ЕЕ 
| 4.60 | 10,79 | 38,64 0,404 | 12,58 | 39,99 | 15661 
0,192 
10,16 | 48,12 
13,25 | 45,00 
0,231 


11,04 | 52,25 0,482 
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Лабораторное моделирование 


Сравнительный анализ критериев подобия, 
полученных в результате обработки данных при 
численном моделировании процесса разрушения 
РВС-30000 м? (натурный объект) и при лаборатор- 
ных исследованиях на модели рассматриваемого ре- 
зервуара в масштабе 1:130, представлен в таблице 4. 

По результатам выполненных исследова- 
ний можно сделать следующие основные выводы. 
Во-первых, из данных, приведённых в таблице 4, 
видно равенство чисел Фруда (Fr = Fr, = idem). Cne- 
довательно, скоростные характеристики потока, об- 
разующегося при разрушении модельного резервуа- 
ра, являются подобными потоку, образующемуся при 
разрушении соответствующего натурного резервуара. 
Во-вторых, во всех рассмотренных случаях числа 
Рейнольдса (в рассматриваемом диапазоне изме- 
нения параметров потока на лабораторном стенде) 
более критического значения (Ке, > Re, = 10000), то 
есть модельный поток, как и натурный, находится в 
турбулентном состоянии. В-третьих, во всех рассмот- 
ренных случаях обеспечивается тождественность 
критерия Струхаля, то есть соблюдается постоянство 
протекания анализируемых процессов во времени. 

В-четвёртых, вследствие того, что скорость по- 
тока воды при проведении лабораторных опытов не 
менее чем в 4,5 раза превышала критическое значе- 
ние (0,23 м/с), то можно утверждать о несуществен- 
ном влиянии на процесс сил поверхностного натяже- 
ния, то есть условие волнообразования соблюдается 
и нет необходимости в увеличении масштаба модели 
исследуемого резервуара. 

Таким образом, соблюдение рассмотренных 
критериев подобия и условий моделирования по- 
зволяет сделать общий вывод о том, что изучаемый 
процесс находится в автомодельной области, а раз- 
работанный лабораторный стенд можно использо- 
вать для нахождения оптимальных геометрических 
параметров различных по конфигурации ограждений 
для РВС и РВСЗС, а также для анализа их устойчи- 
вости к воздействию волны прорыва, образующейся 
при разрушении типовых резервуаров номинальным 
объёмом от 700 до 30000 м?. 
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RESEARCH OF BREAKOUT WAVE PARAMETERS 
AT DESTRUCTION OF TANKS WITH A CAPACITY UP TO 30000 М? 
IN EXPERIMENTAL CONDITIONS 


ABSTRACT 


Purpose. The possibility of using the laboratory- 
experimental method of the research to solve the 
problems associated with finding the optimal geometric 
parameters of different configuration enclosures to 
prevent the spillage fire spread at destruction of typical 
tanks with a nominal capacity up to 30000 плз. 

The object of the study is the process of formation 
and propagation of the breakout wave formed at the 
complete tank destruction. The subject of the study is 
the time-varying speed and height of the wave along the 
spreading path necessary to determine the appropriate 
similarity parameters. 

Methods. Numerical simulation of the destruction 
process of a full-scale vertical steel tank with water is 
carried out in the software complex LSDYNA based on 
the solution of hydrodynamics equations systems and the 
deformed solid in the associated spatial formulation. The 
experimental studies are carried out on the first developed 
laboratory stand with a model tank on a scale of 1:130 
to the full-scale tank. During the stand developing 
conditions of using similarity laws and hydraulic modeling 
rules are observed, and during the experiments the 
measuring equipment was applied which was widely used 
in hydraulic laboratories. 
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Findings. The result of the research shows that 
the process studied in laboratory conditions is in the 
self-similar field (the equality of Froude numbers; the 
identity of the Strouhal criterion; Reynolds numbers on 
the model of more critical value; the condition of wave 
generation is observed). 

Research application field. The laboratory stand 
can be used to find the optimal geometric parameters of 
different configuration enclosures for tanks, as well as 
analyze their resistance to the impact of the breakout 
wave formed during destruction of typical tanks with 
a nominal capacity from 700 to 30000 mê. 

Conclusions. The possibility of using the 
laboratory-experimental method of research Юг 
finding the optimal geometric parameters of different 
configuration enclosures to prevent cascade fire 
development in tank farms of oil and oil products storage 
facilities due to the complete destruction of tanks with a 
nominal capacity up to 30000 m°. 


Key words: experimental simulation, 
destruction, breakout wave, similarity parameters. 
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